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SUMMARY 

The design of a class of special purpose computing machines, which compute by 
counting, is systematically developed. The basis of the design philosophy is to limit the 
basic building elements to three fundamental units and to develop the method of synthesis 
such that these three building elements are represented as operational units. In particu- 
lar, the three basic building elements are (1) the binary rate multiplier, which is a means 
of scaling down a pulse stream to some specified fraction, (2) the counter, and (3) the 
anticoincidence circuit, which is a means of separating pulses arriving at the counter 
simultaneously. The computational errors, that is, rounding-off and truncation errors, 
introduced into the machines when these elements are treated as operational units are 
identified. The method of synthesis is explicitly stated, and a wide variety of machines 
obtained directly from this synthesis are presented. Finally, a series of machines is 
presented for interpolation and extrapolation of a function, which is available only as em- 
pirical data. 


INTRODUCTION 

Computers are usually divided into two broad categories, analog and digital. Analog 
computers represent variables as physical quantities. The solution of a problem in an 
analog computer is attained by constraining a physical model of the problem to be solved 
and measuring the variables. The ability to program a wide variety of problems is 
achieved by having functional components available (e. g. , adders, multipliers, integra- 
tors) and interconnecting them by means of a patchboard. The resulting interconnection 
is scaled to match the desired equation. On the other hand, digital computers represent 
variables as discrete quantities. The usual method of solution of a problem in a digital 
computer is attained by sequencing instructions through the fetch- execute cycle of its 



control unit. Another class of digital computers, known as incremental computers, com- 
bines the parallel functional simplicity and speed of analog computers with the ability of 
attaining computational precision, which is not dependent on precision of measurements. 
Such computers attain a speed advantage over conventional general purpose computers by 
transmitting and processing only partial words in a number of parallel arithmetic organs 
rather than the whole words needed by the fetch-execute cycle. Moreover, the digital 
nature of these computers permits the problem solution to be repeated exactly and there- 
fore does not possess the drift characteristic of analog computers. The incremental com- 
puter, which has commanded the most interest in the literature, is the digital differential 
analyzer, that is, DDA. This computer can be viewed as a digital analogy of an analog 
computer. The usual design practice in each of these machines is to permit them to 
solve a large spectrum of problems. When a computer need arises for a special purpose 
application, this versatility provided by general purpose computers is felt as an added 
cost factor. 

A class of incremental techniques, known in the industry as countup- countdown tech- 
niques, has been used in real time control. Data is represented in these techniques by a 
unitary code. For example, the number 28 is represented by 28 pulses. A function may 
be represented by counting the sequence of pulses in a forward-backward counter or con- 
verting them directly into an analog quantity (e. g. , by a stepping motor) for analog pro- 
cessing. Consequently, when a real time application deals with continuous smoothly 
varying functions, countup-countdown techniques offer a simplicity of hardware that is 
economically more attractive than computing systems designed to handle a large spectrum 
of problems. 

The purpose of this study is to investigate countup-countdown techniques with the ob- 
jective of demonstrating that they can, in fact, be used to generate a wide variety of non- 
trivial functions. This will be done by displaying a circuit that will generate each func- 
tion. However, since the techniques on which this study is based are described in the 
literature only in an ad hoc manner (refs. 1 to 5), only specific circuits will be permitted 
as the basic building elements. In particular, the fundamental units to be permitted are 
(1) the binary rate multiplier (BRM), which is a means of scaling down a pulse stream to 
some specified fraction, (2) the counter, and (3) the anticoincidence circuit, which is a 
means of separating pulses arriving at the counter simultaneously. In order to strengthen 

the previous argument the explicit use of whole word 
adders and subtractors will be avoided completely. A 
succinct recapitulation of the purpose of this study is to 
develop and demonstrate systematically the versatility 
of techniques based on counting for solving sophisticated 
and practical special purpose computer design prob- 
lems. 
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(a) NOR. (b) OR. (c) Delay. 


(e) RST flip-flop. (f) AND. 

Figure 1. - Basic logic elements. 
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The method of synthesis herein will be to describe the principal building elements as 
operational units and then proceed by operational techniques to show how to fabricate the 
various machines. In particular, a first-order difference equation can be represented by 
a counter, and approximate integration can be attained by using a counter in cascade with 
a binary rate multiplier. 

The principal design elements can be described in terms of the logic elements shown 
in figure 1. The NOR element shown in figure 1(a) may have various numbers of inputs. 
The parallel lines shown on one of the inputs in figure 1(f) are included to indicate that 
the AND circuit is intended to act as a pulse gate dependent on the level setting of the 
other line. 


BINARY RATE MULTIPLIER 


A binary rate multiplier is a means of sealing down a pulse stream to some speci- 
fied fraction. A logic diagram of a BRM, which is built of standard logic elements, is 
shown in figure 2(a). This circuit is described in detail in several of the references (e. g. , 


Stage 1 2 k n-1 n 
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<b) Timing diagram. 

Figure 2. - Binary rate multiplier. 
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refs. 2 and 5); consequently, a brief description here will suffice. The input pulse stream 
is applied directly to the binary counter whose value is denoted by x n x n ^ . . . XgXj. 

Each flip-flop of the counter is operated as a trigger. For every two input pulses to a 
trigger two output pulses are produced; one pulse when the flip-flop makes a 0 to 1 tran- 
sition called an a pulse and one when the flip-flop makes a 1 to 0 transition called a 
/3 pulse. The /3 pulse is used to trigger the successive stage of the counter. The a pulses 
are gated through AND gates and mixed through a NOR element to produce the desired 
fraction of the input pulses. This simple mixing technique may be used because the 
a pulses from the various stages are separated in time from each other. This timing 
factor is shown in figure 2(b). 

The quantitative relation of a BRM may be expressed as 

Az = yAx (1) 


where y is the binary number represented by the settings of the AND gates; that is, 


y = . y_ 1 y_ 2 . . . y_ n , Ax is the number of input pulses and Az is the number of output 
pulses. If y remains constant over a Ax interval of 2 n pulses, where n is the num- 
ber of stages of the BRM, then this relation is exact. However, if Ax is less than 2 n 


pulses, this multiplicative relation is only approximate, where the difference between the 
actual output and that given by equation (1) is the rounding-off error in the multiplication. 


It can be demonstrated that this approximation can be improved by increasing the number 


TABLE I. - EXAMPLES OF p-SEQUENCES 


Pulse 


p-Sequence 



Two stage 

Three stage 

Four stage 

1 

10 

100 

1000 

2 

01 

010 

0100 

3 

10 

100 

1000 

4 

00 

001 

0010 

5 


100 

1000 

6 


010 

0100 

7 


100 

1000 

8 


000 

0001 

9 



1000 

10 



0100 

11 



1000 

12 



0010 

13 



1000 

14 



0100 

15 



1000 

16 



0000 


of stages and that it varies with the starting con- 
ditions of the BRM counter. 

Because of the approximate nature of equa- 
tion (1) when Ax is less than 2 n pulses, the spe- 
cific output sequence will be calculated in demon- 
strating specific machines. For these calcula- 
tions, the pulse stream shown in figure 2(b) may 
be displayed in mathematical vector form. This 
will be called the p-sequence. Each position of 
the vector in this sequence represents the possible 
output at a particular pulse time from a stage of 
the BRM. The p-sequence for a two-, three-, and 
four- stage BRM is displayed in table I. 

The p-sequences given in the table assume 
that the BRM counter starting value is zero. If 
another starting value is used, its associated 
p-sequence can be easily obtained. Moreover, if 
an interval greater than 2 n pulses is used, the 
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p-sequence can be obtained by repeating the p-sequence given in the table. 

The sequence of output pulses may be calculated by multiplying bit by bit the 
p-sequence with the respective values of the level settings of the AND gates. This pro- 
cess is illustrated by two examples as follows: 


Example A: 


Example B: 


/ 


100 \ 
010 ' 
100 
001 
100 
010 
100 
\ooo / 


( 01010101 ' 
10111011 
11101111 , 


/ 0 
0 
0 
0 
0 
0 
0 
, 0 



/ioioioiN / 1 \ 

I 0100010 = / 1 1 

\0001000/ 1 
1 

I 1 

w 


The first matrix, in each of these examples, is the p-sequence. The next matrix repre- 
sents successive values of the gate settings. When these two matrices are multiplied to- 
gether by conventional matrix multiplication, the actual result of the BRM output is de- 
veloped along the diagonal of the resultant matrix. This restricted form of matrix mul- 
tiplication is notated by an asterisk and by displaying the result as a vector on the right 
side. The first element of the resultant vector represents the output of the BRM after the 
first input pulse to the BRM' is applied; the second element of the resultant vector repre- 
sents the output of the BRM after the second input pulse is applied, etc. 

The expected output value of example A after the first input pulse is given by equa- 
tion (1) and is equal to 3/8 of a pulse, since y = 3/8 and Ax = 1. The expected output 
value after the second input pulse is the expected output value of the first input pulse plus 
the expected output value resulting from the second pulse. Therefore, the expected out- 
put value of example A is 35/8 pulses for the eight input pulses. As shown by actual com- 
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Figure 3. - Signed four-stage bidirectional counter. 


putation, however, the BRM yields zero output pulses. On the other hand, the expected 
value of example B by equation (1) is 21/8 pulses for the seven input pulses. The pre- 
ceding computation yields seven output pulses. Both of these examples are pathological 
cases in the use of the BRM. The approximate nature of equation (1) can ordinarily be 
expected to yield more realistic results. Some of these results are presented in this 
study. 

The method of synthesis to be presented necessitates that the BRM operate on signed 
quantities. In particular, the level setting of the AND gates and the counter input pulses 
must be signed quantities, and the BRM is to yield a signed pulse output. If the output 
pulses are accumulated in a counter, the sign of the pulse will determine the direction of 
counting. If the output pulses are used to drive a stepping motor, the sign of the pulse 
will determine the direction the stepping motor is to turn. Throughout this discussion 
the signs of various quantities are considered available through level logic. Consequently, 
the output sign can be obtained from the input signs by an exclusive OR circuit. 


COUNTER 

The purpose of the counter in the machines that will be considered is twofold: (1) to 
accumulate the pulses arriving at the counter in order to display the total number of 
counts and (2) to set the levels of the BRM's. In the first application, the counting se- 
quence can be any desired sequence for a terminal device. In many real time applica- 
tions the output pulses may not be accumulated directly but are converted to an analog 
quantity for analog processing (e. g. , by a stepping motor). In the second application, 
the counting sequence must be compatible with the BRM. This general requirement can 
be met by the circuit displayed in figure 3. 

A number is represented in this counter by magnitude plus sign. As had been stated 
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TABLE II. - DOWN -COUNTING 


TABLE m. - UP- COUNTING 


TABLE IV. - COUNTER SIGN 


SEQUENCE 


SEQUENCE 


CONTROL 


Input pulse 

-1 

+1 

+111 

>111 

+110 

-110 

+101 

-101 

+100 

-100 

+011 

-on 

+010 

-010 

+001 

-001 

+000 

-000 

+111 - +001 - -001 

-111 - -001 - +001 


Input pulse 

-1 

+1 

-000 

+000 

-001 

+001 

-010 

+010 

-011 

+011 

-100 

+100 

-101 

+101 

-110 

+110 

-111 

+111 


Sign 

input 

pulse 

Sign 

counter 

Line 

activated 

+ 

+ 

Up 

+ 

- 

Down 

- 

+ 

Down 

- 

- 

Up 


earlier, the signs are represented by level logic. The counter counts down in magnitude 
when the input pulse and counter are opposite in sign and counts up in magnitude when the 
counter and input pulse have the same sign. The circuit is designed so that the a pulses 
are used to count down and the 0 pulses used to count up. There are two representations 
of zero, that is, minus zero and positive zero. When the counter value is at +00 ... 0 
and a -1 pulse arrives, the counter is set to -00 . . . 01. This end correction is ac- 
complished in three steps. The normal sequence first changes the counter value to 
+11 ... 1. The magnitude is then corrected in the second step to +0 ... 01. Finally, 
the sign is changed to -00 . . . 01. The sign is changed last so that the 0 pulses genera- 
ted when the magnitude is corrected do not propagate to the successive stages of the 
counter. In a similar manner to that just given, the counter is set to +00 ... 01 when 
the counter value is -00 ... 0 and a +1 pulse arrives. The down-counting sequence for 
a three-stage counter is given in table n. 

The up-counting sequence utilizing the 0 pulses of the flip-flops is given in table HI. 

Since the signs of both the pulse output of the BRM and counter value are to be pro- 
cessed by level logic, the activation of the up-down line is accomplished by an exclusive 
OR circuit. This is obvious from table IV. 


ANTICOINCIDENCE CIRCUIT 

Pulses arriving at a counter simultaneously must first be separated before they are 
entered into the counter. The circuit that accomplishes this task is called an anticoin- 
cidence circuit. Fundamentally, this circuit necessitates storing each pulse as it ar- 
rives. Each stored pulse is then presented to the counter according to a fixed program. 
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Clock pulses 
C 1 C 2 



Clock pulses 

Figured - Anticoincidence circuit. 

The circuit configuration that can accomplish this task for two inputs is shown in fig- 
ure 4. 

The operation of the circuit given in figure 4 is as follows: If a pulse from input 1 
exists, it is stored in flip-flop 1. If a pulse from input 2 exists, it is stored in flip- 
flop 2. It will be noted that these two inputs can arrive simultaneously. Two pulses are 
emitted by the clock that are separated from each other. If flip-flop 1 had been set by 
input 1, it is reset by the clock pulse C^, which in turn generates an output pulse. If 
flip-flop 1 had not been set, no output pulse will appear in the output. Flip-flop 2 is re- 
set and an output is similarly generated by clock pulse Cg. Since clock pulses and 
Cg are separated, the corresponding output pulses are also separated. 

If either input arrives at the same moment as its associated clock pulse, this circuit 
may fail. However, it is usually the case that the inputs are time dependent on the clock, 
and natural time lags prevent them from arriving simultaneously. In any event, it is 
sufficient for this study just to state that this circuit can be expanded to adjust for such 
failures since the manner by which this can best be accomplished depends intimately on 
the particular application. 

Since the sign of a pulse is processed by level logic, the signs need not be stored be- 
fore they are presented to the anticoincidence circuit. When a pulse is presented to the 
counter, however, its sign must also be presented. This may be simply accomplished 
by shifting the sign level to a flip-flop by the separated clock pulses Cj and Cg. This 
circuit is also shown in figure 4 where S's and S's are the sign levels of the pulses and 
their complements, respectively. 

If more than two inputs arrive at the counter simultaneously, a need arises for a 
circuit other than a simple clock to separate the stored pulses. The design of such a 
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TABLE V. - PULSES GENERATED BY SEVERAL COUNTING 


SEQUENCES 


Counting 

sequence 

Pulses 

generated 

Counting 

sequence 

Pulses 

generated 

Counting 

sequence 

Pulses 

generated 

000 

--a 

000 

--a 

000 

— a 

001 

-qi/3 

001 

-Ol- 

001 

-a- 

010 

--a 

on 

-p 

Oil 

a — 

Oil 

a/3/3 

010 

a — 

111 


100 

--a 

110 

— a 

110 

-0- 

101 

-ctfi 

111 

-/3- 

100 

0- 

110 

-- a 

101 

~/3 



111 

m 

100 

0-- 




multiphase clock is outside the scope of this study. Nevertheless, it is instructive to 
digress and indicate some techniques by which this can be accomplished. First, it is 
noted that a simple binary counting sequence such as the leftmost sequence shown in 
table V will serve this purpose. It will be noted, however, that, while this sequence can 
generate more than two steps, the a and /3 pulses from the various flip-flops are not 
separated, so consequently cannot both be used. 

All the sides of the flip-flops could be used if the counting sequence utilized a unit 
distance code. Such a code would guarantee that not more than one flip-flop would change 
state at any step of the counting sequence. Consider the Gray code counting sequence 
given by the middle sequence in table V. This counting sequence can be used to separate 
as many as six inputs arriving at the counter simultaneously; however, this requires the 
counter itself to go through eight steps. An example of a counting sequence that may be 
used to handle six inputs and yet go through only six steps in the counting sequence is 
given by the rightmost sequence in table V. 


SCHEMATIC REPRESENTATION 

The three circuits described in this section are the principle design elements; how- 
ever, in describing the machines promised in this report, these three circuits will be 
represented as operational units. The advantages to be gained by using operational units 
rather than these circuits are twofold. First, the method of synthesis can be more 
clearly presented. Secondly, a considerable hardware reduction can usually be realized 
when the composite machine is considered because these simplifications arise when all 
the features of these basic circuits are not required. A checklist of the features, which 
when removed would simplify the basic circuits, would include (1) sign control of BRM 
may not be required, (2) counter may not be required both to countup and countdown, 
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Ax 


/^\a Ax 

<!/— 


(a) BRM. 


(b) Scalar BRM. 


Az 


-( 0 ) 


Input 

1 

2 




Output 


(c) Counter. 


(d) Anticoincidence circuit. 


Figure 5. - Schematic diagrams of principle design elements. 


(3) two BRM’s may receive the same input 
pulse stream with the result that one BRM 
counter may be used with two sets of AND 
gates, and (4) the level setting of the BRM 
may be constant with the result that a 
scaling circuit (see ref. 5) rather than a 
BRM may be used. These three circuits 
are, however, sufficient as principal de- 
sign elements. 

The three principal design elements 
as operational units are presented in fig- 
ure 5. The BRM is represented by the schematic diagram shown in figure 5(a). The 
value y in this diagram is less than 1 and is obtained from level logic. The quantities 
Ax and Az are the input and output pulse streams, respectively. The input-output rela- 
tion for this diagram is expressed by equation (1). Alternately, the BRM is represented 
by the schematic diagram shown in figure 5(b) when the value of y remains constant. In 
these cases the BRM may be replaced by a scaling circuit in the final design. Figure 5(c) 
presents the schematic diagram of the counter. The quantity Az is an input pulse stream 
and z is the output that may be used in level logic. When the counter is used to set the 
levels of the gated pulse generators of a BRM, a scale reduction of 2 -n is implied by 
the connection. At times this scale reduction will be shown explicitly by the same dia- 
gram shown in figure 5(b). The initial conditions of a counter may be shown explicitly by 
inserting it in the box, that is, z^. If the counter value is considered to be a function of 
iterative steps, the input-output relation may be expressed by the first-order difference 
equation: 


z (k) " z (k-l) + Az (k-1) 


( 2 ) 


The value of z^ 


in equation (2) in terms of the initial condition of the counter is 


k-1 

Z(k) = z ( o) + 5- (i ) (3) 

Figure 5(d) represents the schematic diagram of an anticoincidence circuit. This circuit 
accepts multiple pulse inputs and produces a single pulse output. The design of this cir- 
cuit is such as to permit the input pulses to arrive simultaneously. At times, however, 
it will be convenient to use this schematic diagram for multiple pulse inputs even if the 
pulses are known to be separated. This usage, therefore, should permit a corresponding 
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Ax J 


M * 


z (0) 


z (k) 


V(k) 


(a) Euler's integration. 


dx 

\ 

y dx ■ 


z 

y — 

i/ 

z (0) 



(b) Integration model. 
Figure 6. - Integrator. 


Entering this output directly into the 


simplification in the final design. Finally, the con- 
vention is adopted in each of these schematic dia- 
grams to use an arrowhead on a pulse logic input or 
output line and not to use an arrowhead on a level 
logic input or output line. 

Approximate integration can be attained by 
using a counter in cascade with a BRM. This con- 
figuration is shown in figure 6(a). If the output of 
the BRM is considered as a function of iterative 
steps, the output at iterative step k is y^Ax. 
counter yields the first-order difference equation: 


z ( k ) = Z (k-D + y(k-i) Ax 


( 4 ) 


K z (0) 


represents the initial value of the z counter, equation (4) may be expressed as 


k-1 

Z W =Z (0) + 5 y W aX <5) 

This equation is recognized as Euler's (rectangular) integration. A model of this pro- 
cess, which is convenient for machine synthesis, is presented in figure 6(b). The devia- 
tion of the results given by the model from that given by equation (5) is the truncation 
error. 

SYNTHESIS (DIFFERENTIAL EQUATION) 

The method of synthesis that will be applied in this section is to express the function 
to be generated as the solution to a differential equation. It has been demonstrated with 
analog techniques that a wide variety of functions can be generated by utilizing only inte- 
grating units and adders (e. g. , refs. 6 and 7). For example, with a mechanical differ- 
ential analyzer, the basic units are a ball-disk integrator and a differential. In the syn- 
thesis of countup-countdown machines, the integrator model of figure 6 and the anticoin- 
cidence circuit, which permits the summation of two pulse streams, will serve as these 
operational units. It should be reemphasized that the principal design elements were 
recognized as entities only for purposes of synthesis, and that the fabrication of the ac- 
tual machines may permit circuit simplification, which may result in reductions of the 
hardware requirements. 

The first step in this synthesis is to express the function to be generated as a 
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differential equation such that the highest order derivative is isolated; that is, 

d m y _ /d m - 1 y dy 

1 7 ? ■ * * 9 9 

dx m V 

The independent variable in equation (6) is represented by a clock. The next step in the 
synthesis is to assume that a circuit has been designed to generate the highest order de- 
rivative. Integrators may then be used to reduce successively the order of the derivative 
according to the equation 



= f dx (7) 

dx k_1 J dx k 

A circuit for generating the highest order derivative, whose existence had previously been 
assumed, may now be developed by the constraint defined by the right side of equation (6). 
The preceding process terminates the design of the basic configuration for generating the 
function. 

The schematic diagram of the machine just designed must then be scaled in order to 
(1) match exactly the defining equation and (2) accommodate the range of variables in a 
finite machine. In particular, the counters are used in the machine configuration to han- 
dle magnitudes that have a finite excursion based on the number of stages. Therefore, 
when a bidirectional counter is used, its magnitude must be such that 


| Counter value | s 2 11 - 1 (8) 

Since the level setting of a BRM must be less than 1, when a counter is used for this pur- 
pose its scale will be reduced accordingly; that is, 

[Level setting of BRM] = 2 -n [Counter value] (9) 

Finally, the scale of both sides of the defining equation, that is, equation (6), must be 
the same. 

This procedure may be reduced to a finite number of steps. A synopsis of these 
steps is as follows: 

Step 1. - Isolate the highest order derivative in the differential equation as shown in 
equation (6). 

Step 2. - Assume the highest order derivative has been generated in a counter. 
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Step 3. - Generate each successive lower derivative by using an integrating unit. 

Step 4. - Constrain the independent variable, the function, and its derivatives ac- 
cording to the right side of equation (6) and connect its output directly to the counter rep- 
resenting the highest order derivative. 

Step 5. - Assign arbitrary constants to the independent variable and its highest order 
derivative. 

Step 6. - Write constraint equations at each counter based on the maximum excursion 
and number of stages. 

Step 7. - Write the constraint equation based on the defining equation. 

Step 8. - Calculate scale factors to satisfy the equations of steps 6 and 7. 

The application of this procedure is illustrated in the following sections in the design 
of specific countdown-countup machines. The first two examples will be machines for 
generating the exponential function and for generating the sine-cosine functions. These 
two machines will be illustrated in detail. 

EXPONENTIAL FUNCTION 


The differential equation 


y = y , y(0) = 1 


( 10 ) 


whose solution is 


y = e 




Figure 7. - Exponential function generator. 


(ii) 

will be used to design the circuit for gen- 
erating the exponential function. The 
design solution for this circuit is presen- 
ted in figure 7. The detailed procedure 
for this synthesis is as follows: 

Step 1. - The defining differential 
equation given by equation (10) is in the 
desired form. 

Step 2. - Assume a circuit has been 
designed to generate the highest order de- 
rivative. This is represented by the line 
labeled y’ in figure 7(a). 

Step 3. - The function y is generated 
by integrating y’. 
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Step 4. - Since by equation (10) the assumed highest order derivative y’ is equal 
to y, then y is directly connected to the line y\ This completes the basic circuit shown 
in figure 7(a). 

Step 5. - This is the first step in the design of the scaled circuit shown in figure 7(b). 
The arbitrary constants A and B are assigned as scale factors to the independent var- 
iable and the highest order derivative, respectively. The interpretation of A is "A 
counts per unit of x. ” The interpretation of B is ”B counts per unit of y'. " Note 
that the scale factor of the counter in figure 7(b) is reduced by 2 -n when it is used to 
set the levels of the BRM. If the value of the counter is By' counts, the level setting of 
the BRM is 2" n By\ 

Step 6. - Constraint equations are written for the counter; that is, 

l 2 ' nAB y|ma* s 2” * 1 

Step 7. - The mechanization of the defining equation is justified by the following con- 
straint equation: 


By* = 2" n ABy 

Step 8. - From the equation of step 7 it can be calculated that 

A = 2 11 counts/unit of x 

The calculation of B depends on the maximum excursion of the variable y accord- 
ing to the equation 


I max 


< 2 n - 1 


This essentially completes the schematic design of the circuit for generating e x . 

. „ „ „ Some choices for A and B based on the equation of 

TABLE VI. - SCALE FACTORS ^ 

step 8 are given in table VI. This series of machines has 
and initial conditions of been simulated on a computer in order to illustrate some 

typical results that can be obtained by machines of this type. 
These results are presented in figure 8. The difference 
equation for the exponential machine may be obtained di- 
rectly from figure 9, which is identical to the circuit shown 
in figure 7 but labeled according to Euler's integration. The 
scale factor associated with the counter is "B pulses per 
unit of y. ” When the counter is at iterative step k-1, its 


EXPONENTIAL MACHINE 


n 

A 

B 

Initial counter 
value 

5 

32 

10 

10 

6 

64 

20 

20 

7 

128 

40 

40 

8 

256 

80 

80 
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value is at y 


(k-lV 


During this iteration, 2 11 Ax pulses arrive at the BRM counter, and 


the BRM puts out By/. Ax pulses. These output pulses are added to the counter to 


'(k-1) 

form the counter value for iterative step k. 
be expressed by the difference equation 


Mathematically the value of the counter may 


B y(k) = B y(k-i) + B y(k-i) Ax < 12 > 

If each clock pulse is taken as an iterative step, Ax = 2~ n and equation (12) may be re- 
written as 


y (k ) = (1 + 2~ n ) y (k _ 1} 



(13) 

Solving equation (13) in terms of the initial 
conditions of y (i. e. , y(0) = 1) yields 

y (k ) = (i + 2 n ) (14) 

This solution is shown in figure 8. 

The difference between the differential- 
equation solution and the difference- equation 
solution is the truncation error of the pro- 
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Ax 

(c) Starting value, 3. (d) Starting value, 4. 

Figure 10. - Output of six-stage exponential machine. 


BRM counter 



cess. The difference between the actual output and the difference- equation solution is the 
error due to round off. The results shown in figure 8 illustrate that the round-off error 
is dependent on the number of stages of the BRM. Figure 10 shows this error for the six- 
stage BRM for a number of different starting values of the BRM counter. This figure il- 
lustrates that the round-off error is dependent on the starting value of the BRM counter. 
By these simple means, measurable improvement can usually be attained in the total error. 

Finally, consider a configuration of this machine built out of standard logic elements 
(see ref. 5) for generating the function. This is shown for a four-stage system in figure 11. 
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Figure 12. - General exponential function machine. 


Since the exponential function is a monotonically increasing function, the counter shown in 
this circuit is a simple forward counter. 

The synthesis and analysis of a countup-countdown machine for the generation of the 
general exponential function 


y = y 0 e 


ax 


( 15 ) 


from the differential equation 


y* = “y, y(0) = y 0 (is) 

follows with only minor modification the design of the machine for generating e x . The 
schematic diagram for this machine is given in figure 12(a), and the logical design is 
given in figure 12(b). 


SINE-COSINE GENERATOR 


The differential equation 


y" = -y, y’(0) = o, y(0) = 1 


( 17 ) 
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(a) Basic design. 



(b) Scaled design. 

Figure 13. - Sine-cosine generator. 


is used to design the sine-cosine generator. The basic schematic diagram and the scaled 
schematic diagram are shown in figures 13(a) and (b), respectively, and may be devel- 
oped systematically as follows: 

Step 1. - The defining differential equation given by equation (17) is in the desired 
form. 

Step 2. - Assume that a circuit has been designed to generate the highest order de- 
rivative. This is represented by the line labeled y" in figure 13(a). 

Step 3. - The function y* is generated by integrating y u , and the function y is 
generated by integrating y\ 

Step 4. - Since y enters the differential equation negatively, the pulses arriving at 
the y counter have a sign change with the result that the output of the counter is -y cor- 
responding to y" (see defining equation). 

Step 5. - Arbitrary constants A and B are assigned as scale factors to the inde- 
pendent variable and the highest order derivative, respectively. The interpretation of A 
is "A counts per unit of x" (i. e. , per radian). The interpretation of B is ”B counts 







per unit of y. " If y' counter is set to 0 and y counter 
is set to - 1 (note that this corresponds to the cosine be- 
ing +1), y’ will countup to generate the sine and y will 
countdown to generate the cosine. 

Step 6. - Constraint equations are written for each 
counter; that is, 

2" lAB |y'lmax s2 ''- 1 

2- 2 “A 2 B|y| max s2“- 1 

Step 7. - The constraint equation is written to justify the defining equation 

By" = -2 _2n A 2 By 

Step 8. - From the equations in steps 6 and 7, the scale factors can be chosen: 

A=2 “’ B |ylmax <2 "- 1 

Some choices for A and B are given in table VH based on the equation of step 8. 
This series of machines has been simulated on a computer, and the results plotted in fig- 
ure 14. 

The truncation error associated with this circuit may be calculated by solving the 
difference equations associated with this circuit. Calling the values of the sine and co- 
sine counters at the k iterative step "BS^" and "BC^ M , respectively, gives the dif- 
ference equations at these two counters: * 

BS (k) = BS (k-l) + BC (k-l) Ax ( 18 ) 

BC (k) = BC (k- 1) ' ^(k-l) Ax < 19 ) 


TABLE VD. - SCALE FACTORS 
AND INITIAL CONDITIONS FOR 
SINE-COSINE MACHINE 


n 

A 

B 

Counter 

Sine 

Cosine 

3 

8 

4 

0 

-4 

4 

16 

8 

0 

-8 

5 

32 

16 

0 

-16 

6 

64 j 

32 

i 

0 

-32 


^qs. (18) and (19) 


may be expressed more concisely as 


1 

(Ax) 2 



o 

where V is the second backward difference. Therefore, the second derivative in 
this machine is approximated by the second backward difference. 
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(23) 


C( k j = (1 + k2" 2n " 1 ) cos 2" n k 

A quantitative evaluation of this circuit is complicated by the fact that it is used to 
generate two functions. One method that seems especially well suited for testing such a 
circuit is to plot one output function with respect to the other function rather than with 
respect to the independent variable. For the sine-cosine generator this is called the 
"circle test" since the resultant figure for a perfect sine-cosine generator would be a 
circle. Moreover, it is possible to study the errors due to round-off independent of those 
due to truncation by comparing the actual output to equation (20) output. For the sine- 
cosine generator a composite plot of the solution to the difference equation may be simply 
obtained by expressing equations (22) and (23) in polar coordinates, with the result that 

p=l + 2~ n_1 0 (24) 

The difference between this equation and the circle represents the truncation error of the 
process and is seen to increase as the spiral of Archimedes. The results plotted in fig- 
ure 14 are compared with the plot of equation (24) in figure 15 by this method. 

OTHER DIFFERENTIAL EQUATION MACHINES 

The hyperbolic -sine - h yperbolic-cosine machine is based on the differential equation 

y" = +y (25) 



Figure 15. - Sine-cosine circle test. 
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p-Sequence 
1 0 0 0 0 
01000 
1 00 0 0 
00 100 
1 0 00 0 
0 1000 
1 00 00 
000 10 
1 0000 
01000 
1 0 000 
00100 
1 0 00 0 
01000 
1 000 0 
0 000 1 
1 0000 
0 1000 
1 00 00 
00100 
1 00 00 
0 1000 
1 0000 
00 0 1 0 
1 0 000 
01000 
1 0 00 0 
00 100 
1 0 00 0 
01000 
10000 
0000 0 


Hyperbolic sine x 

0000000000000000000000000000011111 
0000000000000001111111111111100000 
0000000111111110000001111111100000 
0001111000011110011110000111100111 
0110011001100110100110011001101001 
Output of hyperbolic sine 

0000000000010000010001010101001010 


Hyperbolic cosine x 

1111111111111111111111111111111111 
0000000000000000000000000000000001 
0000000000000000000000001111111110 
0000000000000000001111110000111110 
0000000000001111110000110011000110 
Output of hyperbolic cosine 

1010101010101011101010101010111010 


Figure 18. - p-Sequence calculations for hyperbolic-sine - hyperbolic- 
cosine machine. 



2 dx 



Figure 20. - Square machine. 
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The schematic diagram for this machine is similar to the sine-cosine generator except 
that all outputs from the BRM’s are added to the counters. The basic circuit and the 
scaled circuit for this machine are shown in figure 16. It will be noted from this diagram 
that the values of A and B are defined by the equation set (26) as follows: 

By” = 2 _2n A 2 By 1 


A = 2 n 


y 


Bly I < 2 11 - 1 
l-'maxl 

B I v T I < 2 11 - 1 
I •'max I 


(26) 


The output for a hyperbolic-sine - hyperbolic -cosine generator is plotted in figure 17 for 
a five-stage system where B is chosen equal to 16. It is instructive to display the 
p-sequence calculation from which these results were obtained. These are shown in fig- 
ure 18. 

A series of other useful machines will be illustrated in this section. In particular, 
if two pulse streams du and dv are given, the product uv may be generated by using 
the equation 


duv = u dv + v du (27) 

The basic design for this product machine is shown in figure 19. 

The machine for generating the square of a function is shown schematically in fig- 
ure 20. This machine will generate the function 

y = x 2 (28) 


and is based on the equation 


y’ = 2x (29) 

The square machine is utilized as a subassembly in the machine for generating the 
reciprocal of a function; that is, 


y = 1/x 


(30) 
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The basic design for this reciprocal machine based on the differential equation 

y’ = -y 2 


(31) 


is shown in figure 21. 

The machine for generating the solution to the second-order d if ferential 

y” + 2£u> n y’ + u> 2 y = 0 (32) 

is shown in figure 22. 

The tangent machine is shown schematically in figure 23. This machine is based on 
the differential equation 


y' = 1 + y 2 


(33) 
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A similarity will be noted between this machine and that of the reciprocal machine. 
The square-root machine is based on the solution of the differential equation 


y' = l/2y (34) 

It will be noted by this equation that it will form a subassembly of the reciprocal machine; 
that is, the differential equation 


dz/dy = - 2z 2 


(35) 


is used in order to form the function 


z = l/(2y) 

The basic design of this machine is shown in figure 24. 


(36) 


SYNTHESIS (DIFFERENCE EQUATION) 

Consider the iterative process of successive substitution in the functional equation 

x (k + D = ^k)] (37 > 

where cp(x) is chosen such that the fixed points of cp(x) (i. e. , the points where 
x^ = <p( x j)) are the roots of f(x) = 0. One simple form of cp(x) might be x - f(x), which 
leads to the iterative process 



Figure 24 - Square -root machine. 
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x (k+l) " x (k) ‘ f t x (k)3 


(38) 


A more general form is x - g(x)f(x) which leads to the iterative process 

x (k+l) = x (k) " S[ x ( k )M x (k)] ( 39 ) 

A restriction on g(x) in this latter form is that it has no zeros that are not zeros of f(x) 
and that the multiplicity of its poles at the zeros of f(x) be less than the multiplicity of 
the zeros of f(x) at these points. With these restrictions, it can be readily seen that the 
iterative equation (eq. (39)) has fixed points at the zeros of f(x) (i. e. , x = x - 0). The 
function g(x) in equation (39) is chosen so that the process converges. 

The basic equation of a counter immediately suggests a method for generating an 
equation of the form of equation (39). This method of synthesis is simply to generate a 
pulse stream equal to g(x)f(x) and feed it into a counter. This procedure may be out- 
lined as follows: 

Step 1. - Write the defining equation in the implicit form g(x)f(x) = 0. 

Step 2. - Assume a counter with the value of x and generate a pulse stream equal 
to g(x)f(x)p, where g(x)f(x) is the level setting of a BRM and p is the input to the BRM 
counter for converting this level setting into a pulse stream. 

Step 3. - Feed back the pulse stream generated in step 2 into the x counter. 

Step 4. - Assign an arbitrary constant to each variable represented in the machine. 

Step 5. - Write constraint equations and calculate the scale factors such that these 
equations are satisfied. 


DIVIDE ALGORITHM 

The machine for generating x such that 

x = a/b (40) 


may be designed as follows: 

Step 1. - One way in which equation (40) may be rewritten to put it into the desired 
form is 


xb - a = 0 


( 41 ) 


Step 2. - By assuming a counter value representing x, a pulse stream equal to 
xb - a may be generated (see fig. 25(a)). 
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(b) Scaled design. 

Figure 25. - Divide algorithm. 


Step 3. - The pulse stream generated in step 2 is fed into the counter representing x. 
Step 4. - The schematic of figure 25(a) is redrawn in figure 25(b), and each variable 
is assigned an arbitrary constant. 

Step 5. - The constraint equations may be written directly from figure 25(b), in which 
an iteration is taken at every clock pulse: 



(42) 


(43) 

f 


A = B = C = 2 11 


(44) 
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Equation (44) implies that 


x max — 1 ” 2 


-n 


a <1-2 
max 


-n 


-n 


b <1-2 
max 


and equation (43) may be rewritten as 

x (k+l) = A (k) " * UA (k) 
If equation (46) converges to a fixed point x, 


= X/, v - 2“ n bx/, v + 2“ n a 


( 45 ) 


(46) 


x = x - 2” n bx + 2 -n a 


(47) 


If ^ is the iterative truncation error at iterative step k, that is, 

• e (k) = x - x (k) 


(48) 


then from equations (46) and (47) 


■«(k + i) = ( 1 - 2 "' b > s (k) < 49 > 

This may be written in terms of £ ^ as 

^ (k ) = (! " 2_n b) k ^ (0 ) (50) 

This process will converge if 

Um£^ = 0 (51) 

which implies the condition 

1 1 - 2 -n b | < 1 (52) 

for convergence. Therefore, by equation (52) the process is seen to converge. If the 
sign of the outputs of the BRM are reversed, however, such as shown in figure 26, an 
analysis will show that the process will not converge. 
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OTHER DIFFERENCE-EQUATION MACHINES 

The square- root machine, that is, 


x = Va (53) 

may be designed by finding the zeros of the equation 

x^ - a = 0 (54) 

This machine is shown schematically in figure 27. If the scale factor of x and a are 
both taken as 2°, an analysis similar to that of the divide algorithm shows that the itera- 
tive process 


x (k+l) - x (k) “ 



+ 2 -n a 


( 55 ) 


is generated by the machine. The iteration truncation error for this equation may be 
written as 




( 56 ) 


£ (k+1) - | 1 “ 2 n [ x +x (k)l}' f (k) 
where x is the solution. A sufficient condition for this process to converge is that 

1 1 - 2' n [x + x (k) ] | < 1 (57) 


Since 


x max 


< 1 


(58) 


with the scale factor chosen, the process converges; however, had 

x (k + l) = x (k) ' 2 " n [ a " x fk)J < 59 > 

been chosen as the iterative process, the process would not converge. 

A product machine may also be designed by using this method of synthesis. In par- 
ticular, if 


x = ab 


(60) 


the product may be found by the iterative process 

x (k+l) = x (k) " Pt x (k) - ab ] < 61 > 

This machine is shown in figure 28. If the scale factors for x, a, and b are all taken 
as 2°, the machine generates the iterative process 

x (k + l) = x (k) " 2 ' n t x (k) - ab i < 62 ) 
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( 63 ) 


The iterative truncation error for this machine in terms of £ ^ is 

^(k) = f 1 " 2 ^(0) 


SYNTHESIS (REGENERATIVE CIRCUIT) 

Consider the schematic diagram shown in figure 29. The value of K is bounded such 

that 


|K| < 1 - 2 

The output equation for this circuit may be written as 

dz = K(dx + dz) 

dz _ K 
dx 1 - K 


(64) 


(65) 


As K — 1 in equation (65), the ratio dz/dx — however, equation (64) fixes an upper 
bound on this ratio such that 


K 
1 - K 


= 2 n - 1 


( 66 ) 


dx 


-6 


dx + dz 



dz 


Figure 29. - Regenerative circuit. 


therefore, by using this regenerative circuit, the BRM may act as an amplifier. If large 

amplifications are to be considered, however, 
other factors must be taken into account. In par- 
ticular, suppose a dx pulse arrives at the BRM and 
this generates a dz pulse. The dz pulse is de- 
layed by a gated pulse generator and is fed back to 
the BRM. This in turn may generate another dz 
pulse. This process may be continued depending 
on the value of K; however, each time a dz pulse 
is recirculated, the pulse shape deteriorates. 

For example, if leading edge logic is used, the 
rise time of each pulse will be increased until 
the dz pulse is not sharp enough to be utilized. 
Moreover, enough time must be allowed between 
the dx pulses to permit the maximum number of 
dz pulses. The maximum value of K usually 
utilized in these circuits will in general be less 
Figure 30. - Amplification of regeneration circuit. than that permitted by equation (64). 
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This point is illustrated by the plot of equation (65) in figure 30. If -1 < K < 1/2, at 
most, 1 pulse will be fed back with each input pulse. If 1/2 < K < 1, more than 1 pulse 
will be fed back with each input pulse. Therefore, by fixing the upper value of K, the 
maximum number of feedback pulses may be restricted. 

The method of synthesis in this section is similar to that used in the synthesis by dif- 
ferential equation. In this section, however, the differential equation will involve the 
highest order derivative on both sides of the equation; that is, 


^V = fto,. 

dx m \dx m 



(67) 


In general, the equation is written in this form when the highest order derivative cannot 
be isolated. In particular, if the highest order derivative has a nonconstant coefficient, 
it may be written as equation (67) by adding and subtracting a constant from this coeffi- 
cient (see refs. 8 and 9). The design of the circuit based on equation (67) implies the use 
of the regenerative circuit, since the generation of the highest order derivative involves 
itself. 


SQUARE-ROOT MACHINE 

Consider the generation of the square root 


y = V* 


dy 


C-y 


dx 

T~ 


| + (C - y)dyj 


c-y 


(a) Basic design. 


( 68 ) 


> 



(b) Scaled design. 

Figure 31. - Regenerative-circuit square-root machine. 
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from the equation 


y — = 1/2 
dx 

The first step of this technique is to write equation (69) as 

(y - C + C) & = 1/2 
dx 


and then isolate the highest order derivative as shown in the following equation: 



— + (C - y)dy 
. 2 J 


( 69 ) 


(70) 


(71) 


The basic and scaled circuit for generating equation (71) is shown in figures 31(a) and (b), 
respectively. Following the same technique used to derive a machine based on a differen- 
tial equation gives the constraint equations: 

A dy = B dx + 2 -n A 2 (C - y)dy (72) 


(defining equation for the machine) 


(the counter limiting equation), 
puted by 


A^-ylmax 32 "- 1 < 73 > 

From these equations, the scale factors may be com- 


A = 2“/C 
B=2"- 1 /C 2 




I 


(74) 


where C is chosen such that it satisfies the inequality 


1 


y < i - 2- n 

C max 


(75) 



Figure 32. - Scaled diagram for four-stage square-root machine. 
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Pulses 


Based on this design, figure 32 gives the computed scales and circuit for a four-stage 
regenerative-circuit square-root machine. The output of this machine together with that 
of the desired output is shown in figure 33. 

OTHER REGENERATIVE CIRCUIT MACHINES 

The natural-logarithm machine 

y = In x (76) 

may be designed as a regenerative circuit by the equation 

dy = | + (l-^dy (77) 

The circuit for generating this function is shown in 
figure 34. 

The quotient machine 

z = x/y (78) 

may be designed as a regenerative circuit by the 
equation 

dz = — [(y + C)dz + z dy - dx] (79) 

The schematic circuit for this function is shown in 
figure 35. 


dx 
C 

Figure 34. - Regenerative -circuit logarithm machine. 




x 


I I I I I I i I 

0 2 4 6 8 10 12 14 

Pulses 

Figure 33. - Output of regenerative- 
circuit square-root machine. 
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Figure 35. - Regenerative-circuit quotient machine. 



Figure 36. - Genera! polynomial machine. 


PIECEWISE POLYNOMIAL MACHINES 

When a function to be generated is available only as empirical data (e. g. , such as in 
sampled data systems), the design must be based on generating an approximation function. 
Various classes of approximation functions and techniques for obtaining approximations 
have received considerable attention in the literature (e. g. , see ref. 10). A particularly 
convenient form for approximating a continuous function is that of a polynomial. The gen- 
eral polynomial 


f(x) = 


a m x 


m 


a 


m-1 


ml (m - 1)! 




2 2 

+ — x + a.x + a n 
2 ! 


(80) 


may be generated by the circuit shown in figure 36. However, it is usually the case that 
all the data are not immediately available for generating the function over its entire 
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range, or if it is available, the polynomial needed to meet the accuracy requirements is 
of excessively high degree. In these applications the requirements of the problem may be 
met by using a series of relatively low degree polynomials, where each polynomial is used 
to fit data only in a restricted range. Such machines are called piecewise polynomial 
machines. 

Because of their widespread use in applications (see refs. 4, and 11 to 13) and also 
because they can be realized with relatively simple circuits, this section is devoted ex- 
clusively to describing a series of machines for generating piecewise polynomials arising 
from finite-difference techniques. These machines are grouped into two broad categories 
based on applications, that is, interpolation or extrapolation. Each machine of this ser- 
ies will generate a low order polynomial fitted to data available at equal intervals of the 
argument. In passing from one segment into the next new data are introduced. The form 
of the data in each case is simply generated from the empirical data. 

In order to facilitate the description of the machines in the next two sections, ordinary 
difference notation will be used. In particular, « n is the value of the independent varia- 
ble, where the value of the function is obtained; that is, 

- £ n < 81 > 

2 in 

The quantity 6o> represents the spacing of the independent variable, and A n , A", • . . ,A n 
are the successive differences that may be obtained from lower order differences as fol- 
lows: 


A n f n+l “ f n 


A^ = A . - A 
n n+1 n 




> 


(82) 




1 

1 



J 


In formulating the approximation formula that passes through the given points, it is 
convenient to display these various differences in tabular form as shown in figure 37(a). 
Front this difference table, alternate forms of the approximation formula may be derived 
dependent on the differences that are utilized, that is, on the path through the difference 
table. The paths of some of these formulas are shown in figure 37(b). The form of the 
formula that will be utilized in the following discussion will be such that, in each case, 
the value of the variable x will vary from 0 to 1 in the interval of interest. 
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Gregory-Newton backward formula 



Newton-Sterling formula 
Newton-Bessel formula 


Gregory-Newton forward formula 


(a) Direct-path difference formulas. 


A-i 


A?2 



Newton -Gauss backward formula 


Newton-Gauss forward formula 


(b) Broken-path difference formulas. 

Figure 37. - Difference table and paths of difference formulas. 


p 2 


f(x) 


X 

Figure 38. - Scheme for piecewise 
linear interpolation. 



INTERPOLATION 

A piecewise linear interpolator can be obtained by passing a linear polynomial 
through successive pairs of points of the function to be generated. This scheme is illus- 
trated in figure 38. A linear polynomial is generated that passes through the points Pq 
and Pj. At the point Pj, the point P 2 is added to the scheme and a linear polynomial is 
generated that passes through the points P^ and This procedure may be expressed 

in terms of ordinary differences by the Gregory-Newton interpolation formula; that is, 

f(w n +x6w) =f n + A n x (83) 

The first derivative of this formula is 

Swf^u^ + x6a>) = A n (84) 

This interpolation formula may be generated by the linear polynomial generator shown in 
figure 39(a). 

The corresponding equation and its first derivative for the next interpolation interval 
is given by 


fK +1 + x6co) = f n+1+ A n+1 x 


(85) 
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(b) Function values input data. 

Figure 39. - Machines for piecewise linear interpolation. 


6u,f, K+i + x6w ) = A n +i ( 86 ) 

At the end of the first interval the values of the function and its derivative given by 
equations (83) and (84) are 

f K + M = f n + A n = f n+l ( 87 ) 

6cof’(u> n + 6u>) = A n (88) 

The corresponding values of these two quantities needed at the start of the next interval 
are given by equations (85) and (86) and are 

f< "n + l) * f n+l < 89 > 

6<of'(a> n+1 ) = A n+1 (90) 

By direct computation it may be verified that, in order to proceed from one interval to 
the next, the quantity A^ (which is A n+1 - A fl ) needs to be added to the setting of the 
BRM, and the output (i. e. , the end point of the interpolation interval) need not be modi- 
fied. Since adders have been excluded as basic design elements, however, the same re- 
sult may be attained by transferring A n+1 as the setting of the BRM (since A n+1 = 

A n + A n ). Consequently, the circuit shown in figure 39(a) may be used for piecewise 
linear interpolation of a function by transferring successive first difference as settings 
for the BRM in order to proceed from one interval to the next. 
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The circuit just derived is well suited for an applica- 
tion in which an incremental encoder is used to generate 
the input data. If an absolute encoder is used to generate 
the primary data, the preceding circuit may be adapted for 
this input by using the defining equation for first differ- 
ences, that is, equation (82). This circuit is shown in 
figure 39(b). As in the previous case, only one new piece 
of information must be transferred into the circuit in order 
to proceed from one interpolation interval to the next. In 
this case, however, before the new information, that is, f n+ g, i- s transferred as the set- 
ting of the lower BRM, the value of the lower BRM, that is, f n +1 ’ must be transferred to 
the upper one. 

Two piecewise quadratic interpolators will be derived. The scheme for the first one, 
which will be called the back interval quadratic interpolator, is illustrated in figure 40. 

A quadratic polynomial is generated that passes through points Pq, Pp and Pg. This 
polynomial is used to generate the curve between points Pq and Pp The point Pg is 
then added to the scheme and a quadratic polynomial is derived that passes through the 
points Pp Pg, and Pg. This polynomial is then used to generate the curve between Pj 
and Pg. 

This procedure may also be conveniently expressed by the Gregory- Newton quadratic 


interpolation formula; that is, 

f(u> n + x6u>) = f n + A n x + x ( x ~ ^ A 2 (91) 

The successive derivatives for this formula are 

6u>f’(w n + x6u>) = ^A n - | A 2 j + A^x (92) 

(6w) 2 f"(w n + x6u>) = A 2 (93) 

The corresponding equation and its derivatives for the next interpolation interval are 

f(w n+1 + x6u>) = f n+1 + (i n+1 - i a2 +1 )x < 94 > 

5wf ’K + l + x6 “> - A n + 1 - 1 A n+1 + A n+l x < 95 > 

(6w) 2 f ,, (w n+1 + x6u>) = A 2 +1 (96) 



X 


Figure 40. - Scheme for piecewise qua- 
dratic interpolation (back interval). 
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Consequently, the corrections to be added to the second derivative, the first deriva- 
tive, and the function at the end of the interval in order to proceed to the next interval are 

O 1 O 

A^, - — A^, and 0, respectively. From this formulation, however, addition is required 
in order to proceed from one interval into the next. A formulation of this process that 
leads to the elimination of the explicit adder is to splinter the polynomial given by equa- 
tion (91) into the following two polynomials: 


f (go + x6a>) - f + A x + —x 2 
a v n ' n n g 


(97) 


and 


f b (cu n+ x6o>) = --f x 


(98) 


where 


f(co„ + x6a>) = f (w_ + x6a>) + f,(a> + x6a>) 

II a II D U 


(99) 


The first and second derivatives of equation (97) are 


(6co)f^(a> n + x6w) = A n + A 2 * 


( 100 ) 


and 


(6u)) 2 f”(u> n +x6co) = A 2 


( 101 ) 


The corresponding splintering of equation (94) yields 


f aK + l + x6w > = f n+l + A n + 1 x + x2 


(6co) f 4(a> n+ i + x6cv) = A n+1 + A n+1 x 
(6w) 2 fy(w n+1 + xfiw) = A 2 +1 


( 102 ) 

(103) 

(104) 
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(b) First difference input data. 



(c) Function values input data. 

Figure 41. - Machines for piecewise quadratic interpolation (back interval). 


Consequently, no correction needs to be added to the first derivative in generating 
the function f . Equations (97) and (98) may then be used to design the circuit shown in 
figure 41(a). It will be noted that, in order to proceed from one interval to the next, only 
new second difference data need be transferred to set the levels of the leftmost BRM. 

Based on the defining equation for second differences, that is, equation (82), the 
machines shown in figures 41(b) and (c) may be obtained directly from the machine shown 
in figure 41(a). In these machines, previous first differences and values of the function 
are transferred directly from the lower BRM's to the upper ones before a new first dif- 
ference and a value of the function, respectively, are transferred into the lower one in 
order to proceed from one interval into the next. 
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The scheme for piecewise front interval quadratic interpolation illustrated in fig- 
ure 42 may be derived by use of the Newton-Gauss interpolation formula given in the fol- 
lowing equation: 


'("„ +x6 “)=' n + x V 1 +?fi Y ii4 n-l < 105 > 

If equation (105) is implemented directly, the first derivative and second derivative 

13 3 

must be corrected by adding g ^ n _ i and A n _^, respectively, to these quantities in order 
to proceed from one interval to the next. The explicit need for an adder may be avoided 
in a manner similar to that used in the previous discussion by splintering equation (105) 
into the following pair of equations: 

f aK + x6w > = f n + x A n-1 + \ ^ A n-1 < 106 ) 


f b( w n + x6 w ) = i A n- l x (107) 

Based on this pair of equations, the circuit shown in figure 43(a) may be derived di- 
rectly. The machines shown in figures 43(b) and (c) are adapted from the circuit shown 
in figure 43(a) by using the definition of the second difference given in equation (82). 

A cubic interpolator may be obtained from the Newton-Gauss interpolation; that is, 


f(u> n + x5u>) = f n + x A n + 2 x (x " 1 ) A n- 1 + 1 x ( x ~ l)( x + x ) A n-l (108) 


In this case, however, the discussion will be limited to using this formula for central in- 
terval interpolation only. This scheme is illustrated in figure 44. The points Pq, P 


2’ 


f(x) 


1’ 

and Pg are used to generate an interpolation formula for interpolating between P^ 

and P 0 . The point P. is then added to the scheme and the 


P 3 r 4 

o 


Figure 44. - Scheme for piecewise 
cubic interpolation (central in- 
terval). 


P 2’ 

and 


Pg, and P^ are used to interpolate between 


2 ' 

points Pj, 

points Pg OUU JT 2 . 

Equation (108) may be applied directly to yield a central 

interval cubic interpolator. In this case, however, the third, 

second, and first derivatives must be corrected by adding 
4 14 

A n l , 0, and g A n _ p respectively, to these quantities in ordei 
to proceed from one interval to the next. 

A configuration may be obtained that conforms with the de- 
sign practice of not using an adder by splintering equation (108) 
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(a) Third difference input data. 


... ^n+2» ^n+1 



(b) Second difference input data. 


Figure 45. - Machines for piecewise cubic interpolation (central interval). 


into the following pair of equations: 


f a<“n + “ <n + x ( A n ’ \ a n-l) + ^ A n-1 + 


(109) 


f b (« n+ x6o;) = - | Vi 


(HO) 


Based on this pair of equations, a circuit may be obtained such that the function, its 
first derivative, and its second derivative need not be changed in order to proceed from 
one interval to the next. This circuit is shown in figure 45(a). The circuits presented 
in figures 45(b) to (d) are modifications of this circuit based on the definition of the third 
difference. 


EXTRAPOLATION 

Extrapolation presents an added problem in that the output of the machine must also 
be corrected in order to proceed from one interval to the next. This is illustrated in 


\ 
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W A n+2 



(c) First difference input data. 



Figured. - Scheme for piece- 
wise linear extrapolation. 


Figure 47. - Machine for piecewise linear extrapolator. 


















figure 46 for linear extrapolation. A linear polynomial through Pq and Pj is used to 
extrapolate the values from P^ to Pg. The point Pg is then added to the scheme, and 
a linear polynomial through P^ and Pg is used to extrapolate the next interval. The 
predicted value P* and the new value Pg can be expected to be different. Consequently, 
the output must be corrected for this new value Pg. In order to avoid putting a jump in 
the output function at this point, the scheme that will be employed is to put in the cor- 
rection linearly over the entire next interval. This scheme (as well as that of quadratic 
extrapolation described next) is closely related to the Porter-Stoneman digital filters 
(see ref. 14) and may be extended accordingly. 

The Gregory -Newton backward finite difference formula may be used to design the 
linear extrapolation machine; that is, 

f(a> n + x6u) = f n + A n _ 1 x (111) 

The corresponding formula for extrapolating the next interval is 

f ( w n+l + x5a) ) = f n +l + A n x ( 112 ) 

If these formulas are applied directly, the circuit must be corrected at the end of the 
interval by adding A* ^ to both the function and its first derivative in order to proceed 
into the next interval. This, however, would cause a jump in the output function. This 
jump may be avoided by putting the correction term in the output in a linear manner over 
the entire next interval. The resulting polynomial has the property that its initial value 
corresponds to the end point of equation (111), and its final value corresponds to the end 
point of equation (112). A polynomial that satisfied these constraints may be written as 
follows: 


f (w n+i + x6 co) = (f n + A n _ 1 ) + (A n + A 2 _ 1 )x (113) 

The second difference in equation (113) may be eliminated by using the defining equa- 
tion given by equation (82). This substitution yields the following equivalent equation: 

f Ki+l + x6ct> ) = ( f n + A n-1^ + ( 2A n “ A n -l* X ( 114 ) 

Equation (114) may be implemented to yield the linear extrapolator shown in figure 47. 

The scheme for quadratic extrapolation is shown in figure 48. The quadratic equa- 
tion through points P Q , P^, and Pg is used to extrapolate the data to the point Pg. The 
point Pg is then added to the scheme and can, in general, be expected to be different 
from Pg. The quadratic equation through the points Pp Pg, and Pg is then used to 


47 
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Figure 48. - Scheme for piecewise qua- 
dratic extrapolation. 


extrapolate to the point p£ In order to avoid putting a 
jump in the output when new information is added to the 
scheme, the correction may be put into the output in a 
linear manner over the entire next interval in the same 
manner as that employed for linear extrapolation. 

The Gregory-Newton backward difference formula 
forms the basis of the quadratic extrapolation. This 
formula may be written as follows: 

l K + x6 "> = £ n + x Vl i „ 2 -2 < ll5 > 


The corresponding formula for the next interval is 


f(o, n+ i + tfu) = f n+1 + X A n + Vl 


(116) 


o o 

If equations (115) and (116) are implemented directly, the quantities A„ 9 , 3 A „/2, 
q n- z n- z 

and A n _ ^ must be added to the output function, its derivative, and its second derivative, 
respectively, in order to proceed from one interval to the other. As was indicated ear- 
lier, the jump in the output function can be avoided by putting in the correction over the 
entire next interval. A polynomial that satisfies these constraints (i. e. , has the end of 
eq. (115) as its initial point and the end of eq. (116) as its final point) may be written as 
follows: 


A 2 \ A 2 

+ x6 "> " f n + A n-1 + A n-2 + ( A „ + * 4 n-2] x + "J 1 ** 


(117) 


Substituting the difference relation given by equation (82) into equation (117) yields 

2 


f <“Wl + x£i ") - f n + A n-1 + A n-2 + ( A n-1 + l A n-l ■ 


(118) 


Equation (118) may be splintered into the two equations 


f a<"n+l + x5 "> = Wl + A l 2 + A n-l x + “y 1 ^ 


V“n + l +x6 “> = (!^-l- A n- 2 ) : 

to yield the circuit shown in figure 49. 


(119) 

( 120 ) 
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Figure 49. - Machine for piecewise quadratic extrapolator. 


The circuits shown in figures 47 and 49 may be readily expanded by use of finite dif- 
ference relations (as was done for interpolators) to yield circuits that accept functional 
values and first differences as the primary source of data. 


CONCLUDING REMARKS 

A handful of circuits has been reported in the literature that were designed to meet 
the needs of special purpose digital computer problems arising from real time applica- 
tions. The organization of these circuits is to utilize simple counting techniques as the 
basis for computing. This results in a simplicity of hardware that make them attractive 
for such special purpose applications. The design of these circuits has been examined in 
this study with the objective of (1) 1 ’explaining 1 1 the circuits and (2) generalizing the de- 
sign philosophy such that new circuits may be admitted with the same organization. In 
order to be specific the principal design elements were limited to three fundamental units. 
The elements are (1) the binary rate multiplier, which is a means of scaling down a pulse 
stream to some specified fraction, (2) the counter, and (3) the anticoincidence circuit, 
which is a means of separating pulses arriving at a counter simultaneously. 

These design elements are represented as operational units that may be used to de- 
scribe the machines. Operational techniques are then used as the method of synthesis. 

In particular, a counter is utilized to represent a first-order difference equation and a 
counter in cascade with a BRM is utilized to represent approximate integration. The com- 
putational error, that is, rounding-off and truncation errors, introduced into the ma- 
chines as a result of treating the principal design elements as operational units are iden- 
tified. Having identified these two sources of errors permitted us to obtain better re- 
sults experimentally by two methods: (1) increasing the number of stages and (2) changing 
the round-off error by changing the starting value of the BRM counter. 

The method of synthesis is presented in three parts: (1) expressing the function to be 
generated as a differential equation, (2) expressing it as the fixed point of an iterative 
process, and (3) expressing it in terms of a regenerative circuit that is presented. The 
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method of synthesis is explicitly stated and is satisfactory in that all known circuits may 
be directly obtained from it. A wide variety of other functions is also obtained using 
these synthesis techniques. Many of these examples are illustrated, and in some cases 
actual experimental results were obtained and discussed with the machine. 

A series of machines is presented for interpolation and extrapolation of a function 
that is available only as empirical data. In particular, the function is generated over its 
entire range by a sequence of low order polynomials. Finite-difference techniques are 
used to describe the polynomials. The order of the polynomial is limited to a cubic for 
interpolation, and a quadratic for extrapolation since these seem to be the important 
cases in practice. Nevertheless, these techniques can be easily extended to include 
higher order polynomials. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 21, 1965. 
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